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Abstract. The role of thermal pressure fluctuations in the ejection of tightly packaged DNA from protein 
capsid shells is discussed in a model calculation. At equilibrium before ejection we assume the DNA is 
folded many times into a bundle of parallel segments that forms an equilibrium conformation at minimum 
free energy, which presses tightly against internal capsid walls. Using a canonical ensemble at temperature 
T we calculate internal pressure fluctuations against a slowly moving or static capsid mantle for an elastic 
continuum model of the folded DNA bundle. It is found that fluctuating pressure on the capsid mantle 
from thermal excitation of longitudinal acoustic vibrations in the bundle may have root-mean-square 
values which are several tens of atmospheres for typically small phage dimensions. 
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INTRODUCTION 

It has been emphasized that thermal fluctuations away from equilibrium in small systems assume 
increasingly important roles as the system size becomes ever smaller and approaches molecular 
dimensions [1]. In the case of viral entities, where the phage size may very considerably [2,3], thermal 
pressure fluctuations may become exerted internally by confined DNA that is already tightly packged in 
an equilibrium conformation within a capsid mantle. Such pressure fluctuations can be expected to play a 
signigicant role in the initial packaging and subsequent dynamics of DNA ejection from bacteriophage 
[4,5,6,7,8] when the virus engages a host cell surface receptor and transfers its DNA into the cell. 

Geometrical and topological constraints on the structure and size of viral capsids were discussed in an 
early paper by Caspar and Klug [9] almost fifty years ago. More recent work [2,10,1 1] has discussed 
highly symmetric capsid structure in the form of nanometer-sized protein shells. And viral infectivity has 
been related to DNA length and capsid size [12]. It has also been concluded that very large conformal 
changes may occur in some capsid shells when their viral genomes become tightly packaged [13]. We 
shall adopt a simpler model of viral capsids that contain packaged genomes than has generally been 
discussed [14]. 

CALCULATION OF MODEL ENCAPSIDATED DNA AND CAPSID 

Consider a cylindrical bundle of folded double-stranded DNA bacteriophage segments with circular 
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cross-section diameter L and equal segment lengths L, packaged coaxially and situated within a 
cylindrical capsid cavity of diameter (L+2h) and length L. The bundle of DNA segments is then encased 
by an annular cylindrical shell of thickness h. If the cavity is to have twice the volume [15] of the DNA 
bundle, then geometric considerations lead to the relation h=L / [2(1 +V2)] = 0.2071/L, so that (L+2h) = 
1.414 L is the capsid diameter. Capsids are nanometer sized in general magnitude [10], but vary greatly 
[2,1 1] depending specifically on the DNA genomes they carry within. We shall adopt a representative 
average capsid diameter of 48 nm [5,21] in our calculation of fluctuation pressures exerted between a 
vibrating DNA bundle and the capsid shell to which it is tethered. 

Our discussion of the role of fluctuations rests on a picture of a tightly packaged encapsidated DNA 
bundle that acquires its thermal equilibrium structure through variational minimization of the free energy 
[7] of the capsid-DNA system. Fluctuations of the DNA bundle about its equilibrium conformation takes 
place in the form of thermally excited longitudinal acoustic vibrations. These vibratory DNA excitations 
will be regarded as taking place bounded by capsids which themselves continue to maintain configurations 
close to their equilibrium geometry, and do not participate in thermal excitations on the same more rapid 
time scales as the DNA bundles. The bundles of DNA segments are assumed to be connected internally at 
both ends to their capsid mantles. 

In order to calculate the thermal fluctuations in pressure at the locations where DNA bundles are 
attached to their capsid shells we shall approximate a bundle of folded DNA segments by a continuum 
elastic rod [16, 17]. Longitudinal acoustic vibrations of such a rod fixed at both ends to a static capsid 
inner surface depend on the longitudinal sound velocity within the rod, which in turn depends on the rod's 
macroscopic mass density and elastic properties. These two quantities enter the speed of longitudinal 
sound according to v = V (Y/p) , where Y is the bulk modulus and p is the mass density of the DNA 
bundle. 

We relate these quantities to the microscopic force constant and total genome mass of single double 
-stranded DNA molecules which have become folded into bundles of length L and circular cross-section 
area 7i(L/2) . The elastic response of individual double-stranded DNA molecules has been measured 
[18]. And elastic constants have been used to calculate phase velocities of sound waves [19]. Brillouin 
scattering has been used to determine the longitudinal velocity of sound in B-DNA fibers over a quarter 
century ago [20] with the result v = 1.9 km/s, which we shall adopt for the present calculation. The 
longitudinal sound velocity in a bundle of fibers regarded as a continuum coincides with the sound 
velocity of a single fiber since Y and p have the same ratio. 

We introduce a coordinate system with its origin on the cylindrical bundle's axis where it connects 
with a static capsid mantle. The x-axis coincides with the bundle's axis and extends through it to the other 
end of the tightly packaged bundle on the static capsid' s internal surface a distance L away. Then with 
y(x,t) representing a general longitudinal displacement field, for t > in a continuum bundle we have 

y( x >t) = Zn a n (t) sin(mtx / 2L) , 1 . 

which obeys boundary conditions y(0,t) = and y(L,t) = , where n are even integers and a n (t) are normal 
coordinates for longitudinal wave motion in the folded DNA bundle, regarded as a continuum. We take 
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y(x,t) = for t <0. The total hamiltonian H of the wave system is given by the sum of its kinetic and 
potential energy: 

L 

H = Jo dx [sp/2] [(dy/dt) 2 + v 2 (dy/dx) 2 ] , 2. 

where s = jt(L/2) is the cross-sectional area of the packaged DNA bundle, p is the mass density and 
v is the longitudinal sound velocity. 

We substitute Eq.(l) into Eq.(2) and obtain 

H = (spL/4) Xn K'| 2 + ( v 2 ps7i 2 /16L) £n |a n | 2 n 2 . 3. 

The total energy H depends only quadratically on the a n (t) and their time derivatives. If we now 
adopt a canonical ensemble to obtain the thermal average <H> at temperature T, each quadratic 
term in H has its equipartition average value kT/2 for the ensemble. We obtain: 

< |a n | 2 > = [8LkT/ (nVpsrc 2 )] , 4. 

where the brackets denote the thermal average. 

At the closed end pressure antinode, x = 0, the pressure fluctuation against the constraining internal 
wall of the capsid is 

Ap(0,t) = - v 2 p (dy/dx | x = ) 5. 
for a displacement field y(x,t). We now take the ensemble average < |Ap| > and obtain 

< |Ap| 2 > = [v 4 p 2 n /(4L 2 )] Xnm nm < a n a m > . 6. 

But < a n a m > = <| a n | > 8 nm since in thermal equilibrium the normal coordinates a n are uncorrected with 
respect to their time dependence. A single sum results: 

< |Ap| 2 > = [ v 4 p 2 n /(4L 2 )] Xn n 2 < |a n | 2 > . 7. 
Substituting Eq.(4) into Eq.(7) we obtain 

< |Ap| 2 > = [2v 2 pkT/(Ls)] I n 1 . 8. 

There must be a cut-off limit in Eq.(8) reflecting a short wavelength limit on well-defined longitudinal 
sound waves in the DNA packaged bundle, which corresponds to a highest even mode number N in the 
Eq,(8) summation: 

Z„ N 1 = N/2 . 9. 
3. 



The thermal fluctuation noise pressure on the capsid mantle is then given by: 

< |Ap| 2 >= [Nv 2 pkT/(Ls)] . 
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We obtain a lower bound to the fluctuating pressure magnitude for N=2. Taking M=p(Ls) the 
total genome mass within the capsid, we define the root-mean-square fluctuating pressure magnitudes 
P ms = V(<|Ap| 2 >) and obtain 

Prms> V [2v 2 MkT/(Ls) 2 ] . 11. 

Since M is proportional to the genome length, Eq.(l 1) suggests that the fluctuating 

pressure magnitude at temperature T should be proportional to the square root of the genome length 

divided by the capsid volume whenever the latter is a multiple of the packaged genome volume. 

For our model genome package s = jt(L/2) and the folded cylindrical bundle has half the volume 
of the coaxial caspid cavity. We then have (Ls) = 0.785 L 3 and (L+2h) = L[l+1/(1+V2)] as the capsid 
diameter, which we take as 48 nm [21]. Solving for L we obtain L = 33.94 nm, so that the bundle 

17 3 

volume is (Ls) = 3.07 x 10" cm . Using the virus (b29 as an example [21] we take the genome length 
to be 6.6 - jam . With a base pair separation of 3.4 A we have the number of base pairs in the genome 
equal to 1 .94 x 10 4 for (b29. Multiplying the number of base pairs by the average mass of a base pair, 
1.096 x 10" 21 grams, we obtain the genome mass M = 2.13 x 10" 17 grams. For T = 300 Kelvin and 
v = 1 .9 km/s Eq.(l 1) then gives Prms > 8.11 atm. A virus genome that is sixty times longer, such as 
Mimivirus, would develop an rms pressure of Prms > 63.3 atm. The dependence of ejection forces on 
genome length has been noted in bacteriophage lambda [22] 

In consequence of the above calculations we conjecture that time-dependent fluctuation pressures 
leading to Prms act like a trigger that assists the internal forces exerted on tightly packaged DNA to eject 
the genome. 

I wish to thank Professor Lisa Lapidus for several stimulating conversations. 
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